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17 Abstract 

18 Water ices in surface crusts of Europa, Enceladus, Saturn’s main rings, and Kuiper Belt 

19 Objects can become heavily oxidized from radiolytic chemical alteration of near-surface 

20 water ice by space environment irradiation. Oxidant accumulations and gas production are 

21 manifested in part through observed H 2 O 2 on Europa. tentatively also on Enceladus, and 

22 found elsewhere in gaseous or condensed phases at moons and rings of Jupiter and Saturn. 

23 On subsequent chemical contact in sub-surface environments with significant concentrations 

24 of primordially abundant reductants such as NH 3 and CH 4 , oxidants of radiolytic origin can 

25 react exothermically to power gas-driven cryovolcanism. The gas-piston effect enormously 

26 amplifies the mass flow output in the case of gas formation at basal thermal margins of 

27 incompressible fluid reservoirs. Surface irradiation, H 2 O 2 production, NH 3 oxidation, and 

28 resultant heat, gas, and gas-driven mass flow rates are computed in the fluid reservoir case for 

29 selected bodies. At Enceladus the oxidant power inputs are comparable to limits on non- 

30 thermal kinetic power for the south polar plumes. Total heat output and plume gas abundance 

31 may be accounted for at Enceladus if plume activity is cyclic in high and low “Old Faithful” 

32 phases, so that oxidants can accumulate during low activity phases. Interior upwelling of 

33 primordially abundant NH 3 and CH 4 hydrates is assumed to resupply the reductant fuels. 

34 Much lower irradiation fluxes on Kuiper Belt Objects require correspondingly larger times 

35 for accumulation of oxidants to produce comparable resurfacing, but brightness and surface 

36 composition of some objects suggest that such activity may be ongoing. 

37 

38 


39 



3 


39 Introduction 

40 Explosive volcanism on Earth originates from gas-driven (mostly H 2 O vapor and CO 2 ) 

41 ejection at high temperature of an incompressible fluid from a subsurface magma reservoir to 

42 the planetary surface. High temperatures and plume gas output from volcanic hot spots at Io 

43 suggest a similar propulsive process. Comet outgassing from solar heating on approach to the 

44 Sun is well known, presumably involving release of trapped gases from cometary ices. 

45 Outgassing of N 2 has previously been detected [ Soderblom et al., 1990] by Voyager 2 at 

46 Triton, the largest moon of Neptune, and appears to be similarly driven by solar heating 

47 [Brown et al., 1990; Kirk et al., 1990], but the presence of subsurface fluids is unknown. 

48 Charon now shows potential surface evidence in ammonia hydrates for cryovolcanic flows 

49 from the subsurface [ Cook et al, 2007]. Absence of any ongoing surface change on Europa 

50 [ Phillips et al., 2000] suggests lack in the surface-accessible ice crust environment of 

51 cryovolcanic energy sources. A thick ice crust would likely preclude direct access to 

52 Europa’ s surface of dissolved gases from a putative subsurface ocean indicated by surface 

53 geologic features and induced magnetic fields [Pappalardo et al., 1999]. 

54 Dramatically visible ice plume activity at Enceladus [Porco et al., 2006] now provides 

55 unique opportunities for investigation with repeated observations during the continuing 

56 Cassini mission at Saturn. The measured presence of H 2 O, CH 4 , simple hydrocarbons, and a 

57 yet-unresolved mixture of CO and N 2 [Waite et al., 2006] in the plume gas may indicate 

58 cryovolcanism involving the liquid water reservoir inferred by Porco et al. [2006]. However, 

59 a comet-like process involving only heat-driven release of trapped gases from ice clathrates 

60 has also been suggested [Kieffer et al., 2007]. Although NH 3 is not highly abundant in either 

61 the plume gas [Waite et al., 2006] or in the surface ice [Brown et al., 2006], high temperature 

62 chemistry at the deep core-mantle boundary has been proposed by Matson et al. [2006] as the 
gas and hydrocarbon source. 


63 
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64 Here we propose an alternative process for production of volatile gases that could drive 

65 cryovolcanism on Enceladus and other icy bodies with irradiated near-surface water ices. 

66 This process arises from interaction of continuously produced radiolytic oxidants and 

67 primordial chemistry. First, the unique aspect of this model is that chemical energy for 

68 oxidation arises in the outer ice crusts of planetary bodies from the continuous irradiation of 

69 near-surface water ice by energetic charged particles from the magnetospheric, heliospheric, 

70 and local interstellar space environments. Second, chemical energy and cryovolcanic driver 

71 gases are released by oxidation of primordially abundant reduced “fuel” compounds 

72 including ammonia, methane, and other hydrocarbons that originally accreted with water ice 

73 to form the underlying low density ice mantles of these bodies. Gravitational tides, internal 

74 tectonic activity, and radioisotope decay provide internal heating to increase rates of oxidant- 

75 fuel chemical reactions in cryogenic ice environments. The high exothermic energy yields of 

76 these reactions can potentially make the reactions self-sustaining. The subsurface liquid 

77 reservoir advocated for Enceladus by Porco et al. [2006] would provide an ideal high thermal 

78 gradient environment around the ice-liquid margins of the reservoir to sustain such reactions. 

79 Since reaction rates would increase with temperature in the 80 - 273 K range expected at the 

80 margins of a high-temperature water reservoir, exothermic heat production would raise local 

81 temperatures at reaction sites, and the reactants may concentrate in these locations, explosive 

82 results might be expected as observed. 

83 The primary oxidant, O 2 , in Earth’s atmosphere comes from biological photosynthesis, but 

84 surface oxidants on icy surfaces in the outer solar system are naturally produced in exposed 

85 water ice by radiolysis from irradiation by energetic charged particles and solar ultraviolet 

86 photons. Europa, Ganymede, and Saturn’s rings all have atmospheric environments of 

87 molecular oxygen presumably arising from particle or photon irradiation [Hall et al., 1995, 

88 1998; Johnson et al., 2006], while all three icy Galilean moons have condensed phase O 2 in 

89 the surface ice [Spencer et al., 1995; Spencer and Calvin , 2002]. Ozone is a sensitive proxy 

90 for atmospheric oxygen and is detected at Ganymede [Noll et al., 1996], as well as at the 

91 Saturn moons Dione and Rhea [Noll et al., 1997], Presence of oxidants may be indirectly 
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92 indicated by detection of atmospheric CO 2 as a potential oxidation product at Callisto 

93 [ Carlson , 1999] and of surface CO 2 at Europa [. McCord et al., 1998b], Ganymede [Hibbitts et 

94 al., 2003], and Callisto [ Hibbitts et al., 2000, 2002], Further Callisto spectral observations 

95 and associated models [Strobel et al., 2002; Liang et al., 2005] suggest that atmospheric O 2 , 

96 CO 2 , and CO may be even more abundant on Callisto than the CO 2 first reported by Carlson 

97 [1999]. Cassini spectral mapping measurements [ Buratti et al., 2005] show CO 2 in the ice- 

98 rich surface environment of the small outer moon Iapetus at Saturn. The present paper 

99 explores the possibility that N 2 and CO 2 in the Enceladus plume gas [ Waite et al., 2005] 

100 might arise as oxidation products of radiolytic processes. 

101 The radiolytic oxidant of highest known yield, hydrogen peroxide (H 2 O 2 ), is produced in 

102 low-temperature laboratory measurements [. Johnson et al., 2004; Moore and Hudson, 2000; 

103 Gomis et al., 2004a, 2004b; Loeffler et al., 2006a] at yields of G = 0.1 - 0.4 H 2 O 2 molecules 

104 per 100 eV of ionisation energy deposited by incident primary charged particles. The 

105 measured yield depends on existing abundances of other oxidants and is highest in the 

106 presence of O 2 and CO 2 within the irradiated ice. Although present experimental yields for 

107 direct production of O 2 are highly uncertain [Johnson et al., 2004] by orders of magnitude at 

108 G = 10“ 4 to 10 , even slow accumulation over millions of years in volume ice could saturate 

109 the ice crusts of irradiated icy moons and ring bodies, while also boosting the H 2 O 2 yield and 

110 cumulative concentrations. Sublimation of directly irradiated water ice can also concentrate 

111 the H 2 O 2 product [. Loeffler et al., 2006a], potentially producing explosive [Andrews, 1990] 

112 heating and oxidation transients upon exothermic dissociation of the concentrated peroxide. 

113 Hydrogen peroxide has been detected at the maximally high source rates on Europa’s 

1 14 heavily irradiated surface with surface concentrations ~ 0.1% relative to water ice [ Carlson et 

115 al., 1999] and is possibly present as a UV-absorbing component in surface ices of Ganymede 

116 and Callisto [Hendrix et al., 1999a, 1999b], There are preliminary reports of detection via 

117 near-infrared absorption on Enceladus [Newman et al., 2006, 2007]. The relatively short 

118 lifetimes (four days at Jupiter, twenty days at Saturn) for H 2 O 2 dissociation by solar 
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119 ultraviolet photons make surface detection difficult at sensible depths up to a few hundred 

120 micrometers. Products of more penetrating irradiation at millimeter to meter depths [Cooper 

121 et al., 2001] beyond the UV penetration range -0.15 pm [ Carlson et al. , 1999] would, 

122 however, continue to accumulate. Furthermore, burial of radiation products by water frost 

123 deposition from cryovolcanic emissions may be faster than the production rate. Burial by 

124 meteoritic impact gardening may also protect radiolytic products from photolytic destruction 

125 as suggested earlier for Europa [Cooper et al. , 200 1 ] . 

126 Low abundances of gaseous CH 4 in the atmospheres of Earth, and of Mars [Formisano et 

127 al., 2004; Mumma et al., 2004], illustrate that oxidants and reductants cannot long co-exist in 

128 chemically reactive environments where continuous sources must replenish the observed gas. 

129 The presence of abundant frozen CH 4 on Pluto, Eris, and other Kuiper Belt Objects 

130 [Cruikshank et al., 1997; Brown et al., 2005], and of atmospheric CH 4 on Titan [Niemann et 

131 al., 2005], indicates that oxidized water ices are not presently in direct chemical contact with 

132 the sensible surfaces and atmospheres of these bodies. However, the CH 4 surface layers on 

133 icy bodies, and the dominant CH 4 hydrocarbon in Titan’s atmosphere, could be ejecta from 

134 past and even ongoing cryovolcanism, as arising from occasional interactions of otherwise 

1 35 segregated oxidant and reductant concentrations in the near-surface environment. 

136 Fluid environments are thought to be essential for evolution of life but also enable the 

137 dramatic effects of gas-driven volcanism by the combined gas-piston interaction of 

138 subsurface fluids and gases. Liquid water oceans at tens to hundreds of kilometres in depth 

139 may account for induced magnetic field measurements at the icy Galilean moons of Jupiter 

140 [Kivelson et al., 2004], and a south polar water reservoir may account in part for the plume 

141 activity at Enceladus. Other internal heat sources, e.g. radioisotope decay, may account for 

142 modem presence of fluids even in cases of weak tidal heating, as at Callisto [McKinnon, 

143 2006; McKinnon and Barr, 2006]. If so, even the icy dwarf planets of the Kuiper Belt, many 

144 much larger than Enceladus, might have subsurface reservoirs of the requisite incompressible 
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145 fluids at eutectic temperatures, e.g. 173 K for H 2 O-NH 3 mixtures, to enable active 

146 cryovolcanism. 

147 In this paper we quantitatively model radiolytic gas-driven cryovolcanism for illustrative 

148 cases of Enceladus and elsewhere in terms of the proposed sequence of radiolytic chemical 

149 processes from surface irradiation and oxidant production, to exothermic H 2 O 2 dissociation 

150 and fuel oxidation, and finally to gas production giving rise to plume emissions. Underlying 

151 assumptions in the model are vertical mobility of radiolytic products and underlying fuel 

152 compounds in the south polar ice crust, trace abundances of metal catalysts to trigger H 2 O 2 

153 dissociation even at low temperatures, an unlimited supply of fuel compounds presumably 

154 brought up by rheological flows from deep primordial reservoirs, and a subsurface H 2 O fluid 

155 reservoir that underlies the polar cap. The fluid reservoir provides a steep thermal gradient 

156 around its warm ice margins to accelerate exothermic chemical reactions, becomes 

157 pressurized with gaseous oxidation products, and provides the bulk of plume mass flow 

158 through gas-driven propulsion of the incompressible fluid. Using representative flux spectra 

159 for electron or proton irradiation in selected orbital environments of the outer solar system, 

160 we compute the radiolytic model parameters for icy bodies in these environments under the 

161 above common assumptions to show the full potential range of cryovolcanism that could be 

162 driven by radiolytic gases in diverse outer planet environments from Europa and Enceladus to 

163 the Kuiper Belt and the Oort Cloud. In the final discussion we then characterize the strengths 

164 and weaknesses of this model with respect to those of other models previously published and 

1 65 widely reported for Enceladus. 

166 Enceladus 

167 The water ice [Porco et al., 2006], neutral gas [ Waite et al., 2006], and dust [Spahn et al., 

1 68 2006] plumes of Enceladus provide the second known example of active cryovolcanism and 

169 are also produced, like on Triton, in an icy surface heavily irradiated by energetic particles, 

170 mostly electrons, in the orbital environment of a planetary magnetosphere. We suggest that 
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171 the plume activity at Enceladus could at least partly arise from sporadic interactions of 

172 radiolytic oxidants and primordially abundant reductants brought into direct chemical contact 

173 within the near-surface (tens of meters to kilometers) environment by internal rheologic 

174 processes. Volatile gas products of such reactions, and the primary drivers for cryovolcanic 

175 activity, could include CO 2 and CO, commonly found in outgassing from comets, from CH 4 

176 oxidation, and N 2 from NH 3 oxidation. Plume neutral gas measurements by the Cassini Ion 

177 Neutral Mass Spectrometer (INMS) instrument [ Waite et al , 2006] do not distinguish N 2 

178 from CO, but if N 2 output is significant as now indicated by recent modelling of data from 

179 the Cassini Plasma Spectrometer [Smith et al , 2007], active cryovolcanism and associated 

1 80 outgassing on Enceladus today could be a model for early formation of gravitationally bound 

181 N 2 atmospheres on more massive icy bodies such as Titan [Owen, 2000; Niemann et al, 

182 2005], 

183 Laboratory data [Loeffler et al, 2006b] indicate that NH 3 in thin films of water ice reacts 

184 strongly with radiolytic products of direct energetic particle irradiation to form N 2 gas, but 

185 near- infrared surface measurements [Brown et al, 2006] at Enceladus find no evidence for 

186 NH 3 above trace abundances less than one percent relative to H 2 O. This low surface 

187 abundance is consistent with a similar upper limit 0.5 % reported from neutral plume gas 

188 measurements by the Cassini Ion Neutral Mass Spectrometer [Waite et al, 2006]. Any NH 3 

189 brought to the surface of Enceladus by rheological transport, or as redeposited to the surface 

190 with plume ejecta, is then destroyed rapidly by direct irradiation processes. Absence of NH 3 

191 and other reduced compounds above trace levels also means that oxidants can accumulate 

192 near the surface without immediate losses to oxidation of these compounds. This 

193 circumstance motivates our stratigraphic model (Figure 1) in which oxidant layers are 

194 transported downward from the surface to interact with layers of reduced compounds 

195 similarly transported upward from deep primordial reservoirs. 

1 96 Chemical energy introduced by radiolytic oxidant production at icy surfaces provides an 

1 97 energy source for cryovolcanic activity and could be particularly effective in warm ice and 
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198 liquid near-surface environments suggested [Porco et al., 2006; Spencer et al., 2006] for 

199 Enceladus. The young geologic age from crater counts in the south polar region [Porco et al., 

200 2006], and a current model for polar reorientation by diapirism [Nimmo and Pappalardo, 

201 2006], could indicate ongoing material exchange between the irradiated surface, the 

202 subsurface source layers for the visible plumes, and the deeper interior. Despite assumptions 

203 in other Enceladus plume models to the contrary [. Matson et al., 2006], this vertical transport 

204 cycle apparently does not extend to the moon’s ice mantle boundary with the hot rocky core, 

205 since no Europa-like [McCord et al., 1998a, 1998b] mineral salt species from rocky material 

206 are detected at the surface. Lack of Europa-like non-ice material extrusion or ejection to the 

207 surface, relative to water ice deposition [Verbiscer et al., 2007], is likely why Enceladus 

208 maintains its high visual albedo. Similar considerations may apply to internal transport within 

209 bright Kuiper Belt Objects, perhaps with darker objects being so due to lesser internal 

2 1 0 differentiation during primordial formation as compared to the brighter objects. 

211 Sublimation of the H2O ice matrix in the warm [Spencer et al., 2006] (145 K) south polar 

212 region could produce high H2O2 concentrations [Loeffler et al., 2006a] which could migrate 

213 to the subsurface environment in warm ice flows. High solubility of H2O2 in liquid H2O 

214 would greatly increase downward and lateral transport rates in near-surface melt water or 

215 brine flows and within subsurface water reservoirs. Highly exothermic molecular reactions of 

216 the oxidants could occur in liquid phase with a near-surface reservoir of dissolved reductants 

217 (e.g., NH3) in solution with H2O, or in solid-state phase with hydrate forms of reductants as 

218 illustrated in Figure 1. Mixtures of H2O and H2O2 have lower freezing temperatures than 

219 pure H2O and may be energetically favored in early and present cryospheric environments of 

220 Earth [Liang et al., 2006], Mars [Schulze-Makuch and Houtkooper, 2006; Houtkooper and 

221 Schulze-Makuch, 2006], Europa [Chyba and Hand, 2001], and elsewhere for biological 

222 evolution. 

223 The illustration in Figure 1 shows a variety of oxidant and reductant species, but for this 

224 report we focus on sequences starting with the catalytic dissociation of H2O2 to O2 in the 
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225 presence of a thermal gradient on the ice margin of a subsurface fluid reservoir and of a metal 

226 catalyst X. Exothermic reaction sequences can be ignited in H2O-NH3 or H2O-CH4 mixtures 

227 by the hot O2 from H2O2 dissociation. The principal reaction equations and associated 

228 exothermic heat outputs are as follows: 

229 H 2 0 2 + X->0 2 AEff = 98.2 kJ/mol (1 .02 eV/H 2 0 2 ) 

230 4NH 3 + 30 2 -» 2N 2 + 6H 2 0 A H r = 1359 kJ/mol (14. 1 eV/NH 3 ) 

23 1 CH 4 + 20 2 -» C0 2 + 2H 2 0 AH t = 890 kJ/mol (9.2 eY/CH4) 

232 The activation energy needed to initiate these reactions can be provided in-situ by the 

233 initial H2O2 dissociation or by catalytic reactions of radicals such as OH or HO2. The 

234 subsequent oxidation of NH3 and CH4 respectively forms N 2 , CO, or C0 2 gases as drivers for 

235 cryo volcanic activity. The H 2 0 2 dissociation can be catalysed even in low temperature 

236 environments by trace abundances of dissolved transition metal ions (e.g., Ti, V, Fe, Ni, Cu, 

237 Zn, Mo). Silver or platinum catalysis of H 2 0 2 monopropellent at concentrations of 70 - 98 % 

238 in water is used to produce hot (>873 K) steam mixtures of the H 2 0 and 0 2 combustion 

239 products in rocket propulsion systems [Andrews, 1990], OH and H0 2 radicals are 

240 catalytically produced in the Fenton’s Reagent-like transition metal decomposition of H 2 0 2 . 

241 Iron in the form of FeS04 salt is the typical catalyst for Fenton’s Reagent reactions and the 

242 metal is universally abundant as a potential catalytic agent for Enceladus. Since we are 

243 mainly concerned with visible cryovolvanic activity in the near-surface environments of icy 

244 bodies, these low-temperature catalytic reactions may be more directly linked to such activity 

245 than high-temperature [ Matson et al., 2007] (500 - 800 K) reactions, e.g. for endothermic 

246 NH 3 dissociation (2NH 3 N 2 + 3H 2 ), in the deep interior of Enceladus and other icy moons. 

247 Ammonia has not yet been identified at Enceladus from initial Cassini Ion Neutral Mass 

248 Spectrometer (INMS) measurements of the neutral plume gas [Waite et al., 2006], but 

249 nitrogen as an atomic or molecular ion is certainly present in the magnetospheric plasma of 

250 Saturn at the orbit of Enceladus [Smith et al., 2005; Sittler et al., 2005; Sittler et al., 2006]. 

25 1 Cassini Plasma Spectrometer (CAPS) measurements for radial distributions of nitrogen ions 

252 are consistent [Smith et al., 2007] with a molecular nitrogen source at Enceladus. Apparent 
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253 low abundance <1% of NH3 in the neutral gas does not rule out presence of ammonia in other 

254 forms such as the ammonium ion NH4 + , highly abundant at neutral pH in H2O-NH3 fluid 

255 mixtures and indistinguishable from H2O at the same mass/charge in plasma ion data. 

256 Alternatively, H2O2 oxidation of NH3 may be rapid in the fluid state and consistent with 

257 absence of these species in the plume gas. As depicted in Figure 1, the primary reaction sites 

258 contributing to plume gas production might occur at lower but more sustainable rates in 

259 hydrate [Kieffer et al., 2006] phases within warm ice in boundary regions surrounding a 

260 central liquid reservoir with the main non-H20 components being gaseous products of 

261 oxidation from these regions. 

262 The advantage of gas-driven fluid flow models for cryovolcanism is that modest molar 

263 production rates for gas can lead to much larger mass ejection rates for incompressible fluids. 

264 The measured [ Waite et al., 2006] volatile gas abundances relative to H2O at Enceladus 

2 

265 correspond to a gas-rich subsurface source environment at a total pressure ~ 10“ bars. 

266 Bubbles of hot, 10 K, gas emerging at this pressure into, and rising within, a water reservoir 

267 would expand until pressure and temperature equalize with respect to the ambient hydrostatic 

268 environment. At the one-bar level, nearly equivalent to one kilometer depth in pure 273 K 

2 

269 water at surface gravitational acceleration g ~ 0.1 m s , the mass flow ratio of displaced 

270 water and N2 gas is 800:1 as used for the plume flow rates in Table 1. This ratio 

271 proportionately increases for declining hydrostatic pressure toward the surface until water 

272 depths less than ten meters \Porco et al., 2006] at which sublimation vapour pressure exceeds 

273 hydrostatic pressure. Thus a liquid near-surface environment has an advantage in our model 

274 of high mass ejection for plume fluid relative to gas production. 

275 Enceladus offers the best available test for cryovolcanic models, since in-situ flyby 

276 measurements are available from the ongoing Cassini mission with more planned in the 

277 future. At minimum flow speeds of 240 m s 1 for gravitational escape, the measured or 

278 inferred quantities are H2O mass loss rates: 3 - 90 kg s _l from direct neutral gas 

279 measurements [Waite et al, 2006], 150 - 350 kg s _1 from far ultraviolet absorption [Hansen 
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280 e t al., 2006], and up to 100 kg s _1 from plasma flow deflection by pickup ion mass loading 

281 [Tokar et al., 2006]. Averaged over the south polar cap region of area 72,220 km 2 southward 

282 of 55°S, these limits respectively correspond to kinetic power, delivered to the subsurface 

283 reservoir fluid in our model by oxidation product gases, of 0.001 - 0.04, 0.06 - 0.14, and 

—2 

284 0.04 mWm . Neutral gas data [Waite et al, 2006] also suggest a narrow density spike of 

285 mass flow -0.015 kg s _1 near Cassini closest approach. Overall gas data allow a large range 

286 of spatial and temporal variability in the flow. 

287 As discussed in the next section, the energy fluxes delivered to the south polar region of 

288 Enceladus by magnetospheric electron and cosmic ray proton irradiation, and the heat fluxes 

289 arising from oxidative reactions, overlap the plume flow power requirements. However, the 

290 model gas production rates from radiolytic oxidation are near, or well below, measured lower 

291 limits on plume emissions for the product gases and suggest that relatively high abundances 

292 of gas have already accumulated in a liquid fluid reservoir after depletion of dissolved NH3. 

293 That is, the current radiolytic oxidant input could maintain the kinetic power of the plume 

294 flows but presently be uncorrelated to the cumulative non-H20 gas composition of the fluid. 

295 This suggests a gas-driven cryovolcanic system that is presently in high output phase but that 

296 will eventually decline in output as the reservoir gas pressure is released until more gas 

297 accumulates as we expect from oxidation of reductant fuels. Since the inflowing chemical 

298 energy from radiolysis is continuous, and the potential fuel reservoir from the moon interior 

299 is large, there could be many cycles of high and low activity in the familiar pattern of the Old 

300 Faithful geyser at Y ellowstone National Park on Earth. 

_2 

301 Another measured quantity from Cassini measurements is the ~80 mWm average 

302 surface power of the total 5.8 GW thermal emission from the south polar cap [Spencer et al., 

303 2006] of Enceladus below 55°S. If this emission level were continuous, our estimates on 

304 steady-state radiolytic inputs, as discussed in the next section, would be far too low to 

305 account for it. However, even the globally averaged (area = 798,650 km ) tidal and 

—2 

306 radiogenic internal heating [Porco et al , 2006], each being about 0.5 GW (0.6 mWm ) and 
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307 comparable to the incident energy flux (Table 1) at the global surface from magnetospheric 

308 electrons, fall far short as well. The peak temperature of the measured thermal emission at 

309 145 K [Spencer et al., 2006] is far too low [e.g., Matson et al., 2007] to ignite the oxidation 

310 sequences producing CO 2 and N 2 from hydrocarbons and ammonia, so any direct link 

311 between the heat, water plume, and gas emissions remains highly speculative. Since heat may 

312 be retained for long intervals within a body of sufficiently low thermal conductivity [Kargel, 

313 2006], the current thermal emission could be from an earlier transient pulse of high internal 

314 heating. Conceivably this heat pulse might have created or expanded the putative liquid 

3 1 5 reservoir. The enhanced thermal environment would not directly create cryovolcanic gases 

316 but would do so indirectly through acceleration of catalytic reactions such as we discuss for 

317 the radiolytic model. 

318 Radiolytic oxidant accumulation in the ice crust over tens of thousands of years or longer 

319 could account for transient heat emission and plume activity even at current levels, e.g. if 

320 large clumps of peroxide-saturated ice had recently come into contact with concentrations of 

32 1 ammonia or other fuels as illustrated in Figure 1 . The relatively distinct layers of oxidants and 

322 frost depicted in this figure would break up into clumpier concentrations during downward 

323 rheological transport and could aggregate into distributions of many smaller clumps and some 

324 large clumps over time. Lower freezing temperature of H202-enriched ice could produce 

325 more concentrated pockets of H 2 O 2 over successive local melting and freezing cycles in the 

326 vicinity of the liquid reservoir. Relatively explosive activity may explain the large increase of 

327 February 2004 in neutral oxygen density observed via ultraviolet emissions [ Esposito et al., 

328 2005] as Cassini approached the Saturn system. The inferred doubling of E-ring mass by 

g 

329 5x10 kg and subsequent dissipation to previous mass level over a two-month period is 

330 equivalent to Enceladus outgassing at 100 kg/s over this same period. 

33 1 Energy Input and Gas Production 
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332 The Cassini Orbiter has now flown several times across the jovicentric orbit of Enceladus, 

333 but measurements of the energetic particle flux environment at the orbit of this moon are also 

334 available from the past Pioneer 11 and Voyager 2 Saturn flybys. Representative flux spectra 

335 for the magnetospheric radiation environment of Enceladus and other icy moons of Saturn 

336 have been incorporated from the Pioneer and Voyager data into the SATRAD model [ Garrett 

337 et al., 2005] developed by NASA’s Jet Propulsion Laboratory. The electron flux spectra from 

338 SATRAD are plotted in Figure 2 for the local magnetospheric environments of Mimas, 

339 Enceladus, Tethys, Dione, and Rhea. The respective integrated energy fluxes above 10 keV 

340 for these moons are 2.5, 0.7, 0.5, 0.4, and 0.3 mW/m , so Enceladus’s radiation environment 

341 is of intermediate intensity in the Saturn magnetosphere. 

342 In comparison, the energy flux at Europa in the Jupiter magnetosphere is forty times 

343 higher [ Cooper et al., 2001] even than at Mimas, and yet there is no evidence thus far for 

344 ongoing cryovolcanism or any other surface activity there [Phillips et al., 2000], High 

345 radiation energy flux at icy moon surfaces is evidently not sufficient to produce visible plume 

346 activity. The present cryovolcanic model further requires radiolytic product access from the 

347 irradiated surface to subsurface deposits of reductant fuels such as ammonia and methane. 

348 Cassini magnetometer data on ion cyclotron waves suggest a weaker ion source at Dione by a 

349 factor of 300 as compared to Enceladus [Leisner et al., 2007]. In the context of the radiolytic 

350 model this difference is either due to correspondingly lower abundance of reductants in the 

351 upper ice crust of Dione, or to absence of a liquid subsurface reservoir in the event that 

352 reductant abundances are similar to those of Enceladus. 

353 For Enceladus the SATRAD spectrum, shown again in Figure 3, provides a useful median 

354 sample of the moon’s radiation environment even over the twenty-eight years since the 

355 Pioneer 11 flyby in 1979. Figure 3 additionally shows current low and high limits for flux 

356 distributions of electrons at Enceladus. A collection of other measured MIMI spectra is also 

357 included from several other traversals of the moon’s orbit around Jupiter. These assorted 

358 measurements are mostly well bounded by the upper and lower limit spectra. 
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359 Integrated input energy, H 2 O 2 production, oxidative heat, and N 2 gas fluxes from these 

360 distributions are listed in Table 1 for catalyzed exothermic dissociation of H 2 O 2 to O 2 and 

361 subsequent exothermic oxidation of NH 3 to H 2 O and N 2 gas. For these estimates we use the 

362 maximum measured radiolytic yield, G = 0.4, a reasonable upper limit from laboratory data 

363 [Moore and Hudson, 2000] for an ice crust potentially saturated with CO 2 and O 2 products of 

364 radiolytic chemistry. The plume flow rates in Table 1 are computed for N 2 gas displacement 

365 of H 2 O fluid at the 1-bar pressure (1-km fluid depth) level. 

366 The lower limits on incident energy flux come from the Cassini’s Enceladus I flyby, and 

367 fluxes were substantially higher a quarter century ago during the Pioneer and Voyager flybys. 

368 Subsequent Cassini data show variable fluxes up to the SATRAD level. As a primary source 

369 of neutral gas for the Saturn magnetosphere, high plume output may suppress, via energy loss 

370 through atomic ionisation, the low-energy seed population of energetic electrons in the 

371 magnetosphere. Depletion of low energy plasma electrons in the inner magnetosphere was 

372 earlier found by Voyager as a direct indicator of electron energy loss in the E-ring and 

373 associated neutral gas environment [Sittler et al., 1981], The electron fluxes in Figure 3 are 

374 measured away from the moon and are representative of globally averaged irradiation at the 

375 surface but would vary there with longitude and latitude location. Magnetospheric magnetic 

376 fields [Dougherty et al., 2006] and plasma flow [Tokar et al., 2006] appear significantly 

377 perturbed around Enceladus, so energy flux at the surface may concentrate at higher values, 

378 e.g. by a factor of two, in the polar regions with reduced values at the equator. Enhancement 

379 of more dense or amorphous H 2 O 2 ice in the tiger stripes region [Newman et al., 2006] is 

380 consistent with polar focusing of the magnetospheric particle fluxes. 

381 Comparative incident flux distributions are also shown in Figure 3, along with 

382 corresponding energy flux values in Table 1 , for electrons at Europa [Cooper et al., 2001] 

383 and for cosmic ray proton irradiation of Saturn’s main rings [Johnson et al., 2006] and of 

384 other icy objects [Cooper et al., 2003; Cooper et al., 2006] in the outer heliosphere (40 AU, 

385 heliosheath) and local interstellar medium. Cosmic ray interactions with Saturn’s ring 
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386 material partly accounts for the observed oxygen atmosphere of these rings through total 

387 oxidant production. The dominant oxidant energy source in the rings is presently thought to 

388 be solar ultraviolet irradiation [Johnson et al . , 2006] , and, based on earlier work for Europa 

389 [Cooper et al . , 2001], we estimate - 0.014 mWm as the contribution of this radiation 

390 component at the south polar cap (south of 55°S) of Enceladus. Extrapolation of the cosmic 

391 ray energy flux through Saturn’s planetary magnetic field to Enceladus gives an approximate 

392 surface contribution there well below the contribution limits from magnetospheric electrons 

393 but extending more deeply into the surface. Total UV and cosmic ray energy flux 

—2 

394 contributions at the Enceladus polar cap are then of order 0.02 - 0.03 mWm , one to two 

395 orders of magnitude below the magnetospheric electron input. 

396 The depth profiles of time scales for accumulation of radiolytic dosages computationally 

397 sufficient to convert 90% of all irradiated H 2 O ice molecules to H 2 O 2 at G = 0.4 are shown in 

398 Figure 4 for the electron and proton flux spectra in Figure 3. Within the visibly sensible layer, 

2 

399 less than 1 mm in depth for inferred ice grain sizes of order 10 pm [ Brown et al., 2006], the 

400 time scales are 10 years for Enceladus, 10 years for Europa, and 10 - 4x10 years for the 

40 1 outer heliosphere and beyond. Complete conversion to H 2 O 2 in the visible layer can therefore 

402 occur on moon and minor planet icy bodies, photolytic destruction and other losses being 

403 ignored, over time scales from those of terrestrial ice ages and glacial flow times to the solar 

404 system age. 

405 On Enceladus the 10 -year turnover of the meter-thick icy regolith layer, as approximated 

406 from a Europa model [Cooper et al., 2001] of impact gardening, produces a peroxide fraction 

407 to water ice of nearly 10% in this layer. Paucity of large impact craters in the south polar 

408 region [Porco et al., 2006] indicates continual resurfacing to kilometers in depth. Over the 

409 solar system age, up to ~50 m of water ice at G = 0.4 is then converted to H 2 O 2 on Enceladus 

410 as compared to ~5 km within Europa’s surface and to 20 cm for Kuiper Belt Objects and 

41 1 Saturn’s rings. These values suggest that the upper ice crusts of both Enceladus and Europa, 

412 and much of the sensible and deeper impact regolith layers on KBOs, could be saturated with 
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413 H2O2 and other oxidants. Since the condensed thickness of Pluto’s ~1 0-microbar N2 

414 atmosphere [. Elliot et al., 2003] is equivalent in mass to ~0.1 cm of H2O2, the much thicker 

415 radiolytic accumulation of H2O2 there may have a significant chemical impact on the 

416 atmospheric evolution of Pluto. Higher surface thicknesses ~ 10 cm to 100 m for N2 ice on 

417 Triton [Cruikshcmk et al., 1984], the other site of directly observed cryovolcanic activity, 

418 could easily arise from radiolytic chemistry, since that moon orbits Neptune within a hot 

9 

419 plasma magnetospheric environment. Neptunian energetic electrons and protons provide 10 

420 W energy input to Triton auroral excitation [Krimigis et al., 1990] as compared to our upper 

421 limit of 1.5x10 W from Saturnian electrons at the south polar cap of Enceladus. Greater 

422 distance from the Sun at Triton would further increase relative importance of the local 

423 radiolytic energy source as compared to previously discussed solar energy sources for the 

424 Triton N2 geysers. 

425 Finally, limits on N2 gas mass flux up to 0.08 kg s -1 from H2O2 dissociation to O2 and 

426 NH3 oxidation are given in Table 1 for all cases. Incompressible fluid displacement factors 

3 

427 10 or greater from the N2 mass flux at fluid depths less than 1 km, the one-bar pressure 

428 level, are then sufficient to produce observed total plume mass flows 10* - 10 2 kg s '. 

429 However, the Cassini INMS [Waite et al., 2006] measured upper limits ~0.2 - 6 kg s 1 for the 

430 measured number ratio N2:H20 of 0.04: 1 .0 of Enceladus plume gas are well above even the 

43 1 upper limit of steady-state production from the radiolytic model. 

432 Since the energy and heat fluxes in Table 1 are comparable to limits on the plume kinetic 

433 power, the disparity in gas fluxes suggests a lack of direct coupling between radiolytic gas 

434 injection and plume gas ejection. Together with low measured limits from INMS on 

435 abundance of NH3 in the plume gas, this disparity may indicate that the plume gas is 

436 originating from a gas-saturated subsurface environment in which local abundances of NH3 

437 near the liquid reservoir have recently been depleted by cumulative oxidation. Apparent lack 

438 of H2O2 in the plume gas, if we assume for the moment that INMS has any significant 

439 sensitivity to such reactive species, similarly argues that H2O2 and NH 3 react mainly in the 
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440 boundary regions surrounding the putative liquid reservoir and that only the gaseous 

441 oxidation products persist in that reservoir. Future improvements in laboratory calibration of 

442 INMS with respect to reactive oxidants may change this assessment. Unless the plumes are 

443 unique to the present epoch, the depleted NH3, also including that part lost at the surface by 

444 direct irradiation [Loeffler et al., 2006b], must eventually be resupplied by upward transport 

445 of NH3 - rich ice from the deep interior. 

446 Discussion 

447 The radiolytic model depends on the near-surface production of oxidants that should be 

448 present near the surface of icy bodies of the outer solar system, wherever water ice is 

449 irradiated by energetic particles and ultraviolet photons from the magnetospheric, 

450 heliospheric, or interstellar space environments of these bodies. Chemical energy 

45 1 accumulates in these oxidants over a vast range of time scales in different magnetospheric, 

452 and heliospheric, locations. This energy can be released when the oxidants come into contact 

453 with reductant fuels such as NH3 and CH4 in the presence either of high temperatures or 

454 catalytic materials Hydrogen peroxide is produced by ice radiolysis and photolysis in low- 

455 temperature environments and has the advantage of exothermic, potentially explosive, 

456 dissociation to O2 in the presence of iron and other commonly abundant metal catalysts. No 

457 laboratory measurements are yet available on metal-catalyzed dissociation rates at the 

458 observed temperatures of 80 - 145 K on Enceladus but these rates would rapidly increase in 

459 the warmer ice surrounding a liquid water reservoir. This initial exothermic reaction could 

460 then initiate oxidation of reductants to produce volatile gases driving cryovolcanic activity. 

461 Trace hydrocarbon species, such as the INMS-detected [Waite et al., 2006] acetylene and 

462 propane at Enceladus, could be produced in the high temperature sequences initiated by H2O2 

463 dissociation. These species may take part in further reactions, potentially with explosive 

464 results, as found by Benit and Roessler [1993] for proton-irradiated acetylene frost at 77 K, 

465 even at the low end of the Enceladus temperature range. 
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466 Visible manifestations of cryovolcanism include emissions of ice grains, neutral gas, 

467 pickup ions, and dust, but the generally bright reflective surfaces of other icy bodies, as in the 

468 Kuiper Belt, are also highly suggestive of sustained cryovolcanism. New observations of 

469 ammonia hydrate on the surface of Charon in the Pluto system [Cook et al., 2007] suggest 

470 that cryovolcanism is active there. In the Saturn magnetospheric environment the Enceladus 

471 ice grain output sustains the E-ring population which then globally bombards and brightens 

472 [Verbiscer et al., 2007] the surfaces of Enceladus and other moons. Cyclic activity at 

473 Enceladus could produce stratigraphic layers (Figure 1) of high oxidant accumulation during 

474 low activity intervals and of high albedo water frost from direct plume output and surface 

475 bombardment by E-ring grains during high activity phases. Currently high plume and thermal 

476 emission suggests that Enceladus is now in its most active phase, and that any non-ice 

477 materials, including abundant oxidants expected from continuous magnetospheric irradiation, 

478 are now mostly buried from view by water ice frost from plume ejecta. 

479 Enceladus outgassing into Saturn’s magnetosphere, and resultant electron energy loss in E- 

480 ring dust and neutral gas, should lead to anti-correlation of plume output with the fluxes of 

48 1 energetic and lower-energy suprathermal electrons near and beyond the orbit of Enceladus. 

482 The large decrease in low energy (< 0.5 MeV) electron and energetic (> 0.5 MeV) ion fluxes 

483 inward from the orbit of Rhea (8.75 Rs) to that of Enceladus (3.75 R$) was first discovered 

484 by Pioneer 11 measurements [McDonald et al., 1980; Simpson et al., 1980; Trainer et al., 

485 1980; Van Allen et al., 1980], The subsequent Voyager measurements found dropouts in 

486 suprathermal (0.03 - 6 keV) electron [Sittler et al., 1981], ring current ion [Connerney et al., 

487 1981; Acuna et al., 1981], and energetic (> 500 keV) ion [Krimigis et al., 1981, 1983; and as 

488 reinterpreted by Paranicas et al., 2004] fluxes. 

489 These inward flux decreases are now clearly attributable to Enceladus neutral gas and E- 

490 ring ice grain interactions a quarter century after the first suggestion of a potential E-ring 

491 interaction by Thomsen and Van Allen [1979] in a report published just after the Pioneer 1 1 

492 flyby. In comparison, the Pioneer and Voyager energetic electron fluxes increased with no 



20 


493 evident E-ring interaction toward Saturn, until cutoff by the main rings, from source regions 

494 in the Enceladus - Rhea region and beyond. Generally reduced fluxes of energetic electrons 

495 near the orbit of Enceladus, as recently measured by the Cassini MIMI instrument (Figure 3), 

496 suggest that the currently high plume output may be attenuating the supply of low-energy 

497 electrons feeding via magnetospheric acceleration processes into the energetic electron 

498 population. If so, the current rate of surface oxidant production is now low but will again 

499 increase if plume output subsides and the low-energy electron source of the energetic 

500 electrons at Enceladus is no longer attenuated in diffusing inward through the E-ring region. 

501 Negative feedback of Enceladus plume activity to the electron source could produce cyclic 

502 high phases of “Old Faithful” activity initiated after long periods of high oxidant 

503 accumulation rates during low activity phases. Cyclic plume activity elsewhere may be 

504 confined to planetary magnetospheric environments in which the cry o volcanic emissions are 

505 large enough to attenuate source fluxes of energetic particles producing the surface oxidants. 

506 In heliospheric environments the activity may be more episodic and driven by subsurface 

507 circulation of reductant fuels to the oxidant-saturated near-surface water ice environment. In 

508 all cases, interactions of clumps of highly concentrated oxidants and reductants, as in Figure 

509 1 for Enceladus, could produce transiently high plume activity and heat output after slow 

510 accumulations over much longer time intervals, e.g. thousands of years for Enceladus as 

5 1 1 compared to the quarter-century of our present Voyager-Cassini measurements for the E-ring 

512 and now for Enceladus. In comparison, long term storage of heat, e.g. due to low thermal 

513 conductivity in the ice mantle [Kargel, 2006], can be theorized for alternative theories 

514 invoking episodic tidal heating, but we have yet no independent evidence for such episodes in 

515 the surface features and chemical composition of Enceladus. 

516 Enceladus has the highest density, 1.6 g/cc, of all major Saturn satellites excluding Titan, 

517 but the sensible surface does not show composition consistent with continuing chemical 

518 contact to the requisite rocky core [ Porco et al., 2006] underlying the outer ice mantle at this 

519 density. There is no evidence from near-infrared imaging spectroscopy [Brown et al., 2006] 
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520 for surface presence of Europa-like [McCord et al., 1998a, 1998b] sulfate hydrates or other 

521 non-ice species that might directly convect through warm ice crust layers to the surface or 

522 become manifested through the plume outflow. We suggest that this poses a significant 

523 problem for the widely publicized theory of Matson et al. [2007] that the plumes are a 

524 manifestation of high temperature (> 500 K) ammonia chemistry and resultant gas production 

525 deep in the moon interior at the core-mantle boundary. A deep source is also difficult to 

526 connect to plume variability strongly suggested by the large transient increase [Esposito et 

527 al., 2005] in magnetospheric neutral gas during Feb. 2004. Such transients are more likely to 

528 arise from plume source dynamics closer to the surface. 

529 We acknowledge the suggestions of Kargel [2006] and Kieffer et al. [2006] that gas 

530 clathrate decomposition in the presence of thermal gradients could contribute to Enceladus 

531 cryovolcanism, but production of the enclathrated gases then remains to be explained. If 

532 formed during original accretion and differentiation of Enceladus, these gases would have 

533 long since escaped, since this moon is too small to retain an atmosphere like that of Titan 

534 [Owen, 2000], Sequestration of the ammonia in rocky ammonium minerals was suggested as 

535 one alternative [Kargel, 2006], but lack of surface compositional expression for continuous 

536 ice mantle circulation from the rock core boundary provides no current support for this 

537 alternative. The radiolytic model instead postulates that primordial ammonia, methane, and 

538 other hydrocarbon sources of cryovolcanic gases are abundant within the ice mantle, e.g. in 

539 the form of hydrates, and that gas production occurs in the warm ice basal margins of a near- 

540 surface liquid reservoir on contact with oxidants from the irradiated surface. 

541 Even at the relatively high temperature of pure liquid water, the spontaneous dissociation 

542 of H 2 O 2 to O 2 must be catalysed in the radiolytic model by non-ice contaminants, e.g. 

543 universally abundant Fe, originating either from the rocky mineral composition of the deep 

544 interior or from meteoritic bombardment. The absences of detectable minerals on the surface 

545 and of O 2 in the atmospheric environment argue that this dissociation cannot occur very near 

546 the surface, and that radiolytic H 2 O 2 can continuously accumulate and concentrate in the 
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547 near-surface ice environment. If the near-surface environment undergoes cyclic melting and 

548 freezing phases, the non-ice materials will accumulate as brines at lower depths below the 

549 melt layers. These contaminants will similarly have low abundances within a liquid 

550 subsurface reservoir and likely contribute only at trace levels to composition of cryovolcanic 

551 plumes. For the present model we assume that sufficient trace abundances of Fe and other 

552 metals are mixed with continuously upwelling fuels, e.g. hydrates of NH 3 and CH 4 , to trigger 

553 H 2 O 2 dissociation at the basal thermal margins of the fluid reservoir as depicted in Figure 1, 

554 and that these metals do not limit the accumulation of H 2 O 2 above the reservoir. This vertical 

555 gradient in catalyst density allows the requisite inverse gradient in oxidant accumulation to 

556 occur and sets the stage for the exothermic reactions leading to gas production and fluid 

557 ejection to form cryovolcanic plumes. Initial vertical segregation of oxidants and fuels may 

558 be essential to the explosive form of cryvolcanism potentially accounting for the spectacular 

559 plumes of Enceladus. 

560 Saturn kilometric radiation (SKR) provides a new diagnostic for limits on magnetospheric 

561 mass loading by ions originating from Enceladus plume gas and ice grains. During the 1980 - 

562 1981 flybys the two Voyager spacecraft discovered these radio emissions with a periodicity 

563 of 10 hours 39 min 24 ± 7 s, then thought to arise from intrinsic rotation of Saturn [Desch 

564 and Kaiser, 1981]. Later Ulysses radio observations from 1994 to 1997 showed changes in 

565 SKR periodicity up to 1% relative to the Voyager observations. Initial Cassini measurements 

566 determined a more precise value of 10 hours 45 min 45±36 s [ Gurnett et al., 2005] for an 

567 increase of 6.35 minutes (1%) since 1980 - 1981. Slippage of the coupled magnetosphere- 

568 ionosphere system with respect to planetary rotation, due to magnetospheric plasma loading 

569 by Enceladus outgassing and subsequent ionisation of plume neutrals, is proposed [Gurnett et 

570 al., 2007] to account for the radio SKR period. A theory of rotationally-driven plasma 

571 convection [ Goldreich and Farmer, 2006] estimates the mass loading for the 1% SKR period 

572 variation at 10 kg/s radial mass outflow from Enceladus and the E-ring. Although zonal 
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573 variations of Saturn’s atmospheric rotation speeds presently preclude independent 

574 measurement of intrinsic planetary rotation, differential changes in SKR period could be 

575 monitored as a linear measure of Enceladus plume output. Stability of this period within 1% 

576 does imply relative stability over the last quarter century but does not limit longer term 

577 stability. Continued monitoring of SKR periodicity, as compared to Enceladus plume activity 

578 and magnetospheric particle fluxes, is needed to further quantify Enceladus two-way 

579 coupling to the magnetosphere. 

580 The current plume flow limits of 3 - 350 kg/s for Enceladus project over the solar system 

581 age to the cumulative loss of 0.4 - 46 % from the moon’s total mass. For spherical symmetry 

582 a flow rate of 70 kg/s extrapolates to complete loss of the nine percent of total moon mass 

583 southward of 55°S. Enormous redistribution of mass would have globally occurred within the 

584 deep interior if the plumes had been continuously active at such rates. A far lower limit on 

585 cumulative mass loss of 0.03%, corresponding to 0.2 kg/s in flow rate, comes from the 

586 observed ~ 0.5-km depression of the south polar cap region [ Porco et al., 2006], although 

587 much of the plume mass has likely returned to the surface. Thus the average flow rate may 

588 have been even much lower even than 3 kg/s, and there could have been long intervals for 

589 accumulation of stored radiolytic energy to fully account for all aspects of current plume flow 

590 and polar cap heat emission. The present epoch of high activity, manifested in the E-ring 

591 during the Voyager era of 1980 - 1981, and now seen directly by Cassini, may be a short- 

592 lived and periodic event over the history of Enceladus. An apt analogy may be Old Faithful 

593 with variable output eruptions over a few minutes at intervals of about ninety minutes. It is 

594 therefore now advisable to calibrate remote long-term indicators of plume activity, potentially 

595 including SKR emission that be monitored from Earth, to the ongoing Cassini measurements 

596 within the Saturn system. 

597 Oxidant chemistry induced by cosmic ray irradiation could drive resurfacing on the large 

598 bright Kuiper Belt Objects such as Eris and more generally account for the diversity of 

599 surface brightness and color [Jewitt and Luu, 2001; Doressoundiram et al., 2002] among the 



24 


600 icy dwarf planets of the outer solar system. Water ice mantles should be common on KBOs 

601 and deep water-ice absorptions are directly seen on some objects such as the 2003 ELg) 

602 collisional family [Brown et al., 2007; Barkume et al., 2006] and Charon [Brown and Calvin, 

603 2000], Although N 2 and CH 4 ices are dominant on larger bodies such as Pluto [Brown and 

604 Calvin, 2000] and Eris [Brown et al., 2005], surface layers of these volatile molecules would 

605 not be retained on smaller objects [Brown and Calvin, 2000] and, if present, would have to be 

606 renewed by outgassing from the subsurface, e.g. by cryovolvanism. Oxidant production from 

607 cosmic ray irradiation of visible surface or near-surface water ice must therefore be 

608 accounted for as one chemical energy source for potential cryovolcanic activity of these 

609 objects. The respective Voyager 2 and Cassini flyby observations have given the first close 

610 looks at active cryovolcanism on Triton and Enceladus. New Horizons is our next 

611 opportunity to search for such activity at Pluto and Charon during the April 12, 2015 flyby of 

612 this icy dwarf planet system. 

613 Finally, the suggested vertical segregation of near-surface radiolytic oxidants and deeper- 

614 lying reductant fuels on irradiated icy bodies potentially provides chemical energy resources 

615 for future exploration and habitation of the outer solar system and beyond. Since heliocentric 

616 orbits of many of these bodies, such as Centaurs, are dynamically unstable, oxidant and fuel 

617 mining can be envisaged on icy objects having more recently entered the inner solar system 

618 after billions of years of radiolytic processing in the Kuiper Belt or Oort Cloud regions. If 

619 there is ever the need to divert such an object from impact with the Earth, this potentially 

620 accessible chemical redox energy accumulated within the body could become a critical 

621 resource. In these respects, if the radiolytic model proves to be correct, Enceladus may be 

622 viewed not only as an object of intense scientific interest but also as a natural model for 

623 astroengineering. Sagan [1994] suggested utilization of icy bodies as spaceships to other 

624 stars, and Enceladus may be a spectacularly visible model for fulfilment of this vision. 


625 
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Table 1. Radiolytic Cryovolcanism Model Parameters for Enceladus 



Particle 

h 2 o 2 

Chemical 

n 2 

Plume 

Body - Model 

Energy 

Number 

Heat 

Source 

Fluid 

Flux 

Flux* 

Flux 

Flow 

Flow 


mWmf 2 

(cm 2 -s)~' 

mWnrf 2 

kgs -1 

kgs' 1 

Europa - e 

100 

1.7x10" 

5.4 

3.8 

3000 

Enc - high e 

2.1 

3.6xl0 9 

0.11 

0.08 

64 

Enc - mid e 

0.7 

1.2xl0 9 

0.039 

0.027 

22 

Enc - low e 

0.1 

1.7xl0 8 

0.0054 

0.0038 

3.0 

Enc - CRP 

0.002 

3.4xl0 6 

0.00011 

7.7x1 0~ 5 

0.062 

Rings - CRP 

0.0005 

8.5x10 s 

2.7x1 O' 5 

1.9x1 O' 5 

0.015 

40AU SolMin 

0.003 

5.2xl0 6 

0.00016 

0.00012 

0.096 

- CRP 






HS - CRP 

0.0035 

6.0xl0 6 

0.00019 

0.00013 

0.10 

LISM - CRP 

0.0083 

1.4xl0 7 

0.00044 

0.00032 

0.26 


845 *Computed at maximum radiolytic yield G = 0.4 for H 2 O 2 . 

846 Abbeviations: Enc - X - e (high, middle, and low electron flux models for Enceladus), CRP (cosmic ray 

847 proton), Ring (Saturn’s main rings), SolMin (solar minimum), HS (heliosheath), LISM (local interstellar 

848 medium). 
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851 Figure Captions 

852 Figure 1. Illustration of radiolytic gas-driven cryovolvanism for potential subsurface liquid 

853 water reservoir near the surface of Enceladus. A continuous rain of energetic electrons 

854 (yellow arrows) drives radiolysis and saturates the upper ice surface with oxidants (blue), 

855 mostly H2O2 but with mix of other species. Transient plume activity ejects ice grains falling 

856 back to the surface (white ballistic curves) and depositing multiple layers of water frost 

857 (white) interspersed with oxidant concentrations accumulating during lower plume activity. 

858 Ice upwelling continuously supplies fresh “fuels” (red), e.g. hydrates of CH4 and NH4, from 

859 the deep interior. Exposure to increasing temperature within the ice margins of the fluid 

860 reservoir initiates sequences of exothermic reactions from H2O2 dissociation to fuel 

861 oxidation, resultant gas production, and fluid reservoir heating. Percolation of expanding hot 

862 gas bubbles from the margins into the incompressible water fluid becomes the driving force 

863 for upward movement (white arrows) of gas-saturated fluid to form the plumes. Interactions 

864 of lower (left side) or higher (right side) concentrations of oxidants and fuels produce 

865 correspondingly less or more gas. 

866 Figure 2. Representative differential flux spectra of Saturn magnetospheric electrons at the 

867 equatorial jovicentric orbits of the icy moons Mimas, Enceladus, Tethys, Dione, and Rhea as 

868 determined from the SATRAD model [Garrett et al., 2005]. 

869 Figure 3. Differential flux spectra of energetic electrons (red curves) at Enceladus and 

870 Europa, respectively from the planetary magnetospheres of Saturn and Jupiter, and of cosmic 

871 ray protons (black curves) irradiating icy bodies within and beyond the heliosphere. The 

872 thickest solid Enceladus curve is from the SATRAD model. The two medium-thickness 

873 Enceladus curves are respectively upper and lower limits from Voyager and Cassini/MIMI 

874 electron data. The upper limit comes from Voyager measurements compiled by Maurice et 

875 al. [1996], and the lower limit is from Cassini electron flux measurements by the 

876 Magnetospheric Imaging Instrument (MIMI) experiment [Krimigis et al., 2004] near the 
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877 Jupiter-fixed longitude of this moon during the 9 March 2005 flyby. The limiting flux curves 

878 for Enceladus are extrapolated into the SATRAD curve at higher energies. The thinnest red 

879 curves are MIMI measurements for other selected periods of Enceladus orbit crossings: Julian 

880 days 48, 195, and 266 of 2005. A small peak at 0.03 MeV in the spectrum from day 48 is 

881 from a transient injection event. The Europa flux spectrum [Cooper et al., 2001] is a 

882 composite of Galileo Orbiter data below 1 MeV from the Energetic Particle Detector 

883 [ Williams et al., 1992] and Pioneer-Voyager model data [Divine and Garrett, 1983] at higher 

884 energies up to 40 MeV. The Enceladus and Europa model spectra are extrapolated for surface 

885 interaction modelling as power laws to 1000 MeV from the lower energy measurements. 

886 Proton flux spectra are also shown as derived for cosmic ray protons (four dashed curves) in 

887 the outer supersonic heliosphere within the classical Kuiper Belt near 40 AU and in the outer 

888 heliospheric boundary regions [Cooper et al., 2006] of the heliosheath and the local 

889 interstellar medium. The two 40-AU spectra are respectively for minimum [Cooper et al., 

890 2003] and maximum solar activity, the latter data being provided from previously 

89 1 unpublished Voyager data of co-author M. E. Hill. 

892 Figure 4. Surface depth profiles in H 2 0 ice for time in years to accumulate chemically 

893 significant dosages of 60 gigarad (~ 110 eV per irradiated H 2 0 molecule) from total 

894 magnetospheric electron irradiation of Enceladus and Europa as compared to trans-neptunian 

895 objects. The two thin solid curves show the partial contributions of higher energy electrons at 

896 1-10 MeV and 10-100 MeV to the Enceladus total irradiation at 10 keV to 1000 MeV, as 

897 compared to the total irradiation profile for Europa (dot-dot-dash curve). Comparative 

898 profiles (dashed curves) are also shown for cosmic ray proton irradiation of objects near 40 

899 AU near solar minimum [Cooper et al., 2003] and on highly eccentric orbits [Cooper et al., 

900 2006] passing through the heliosheath and into the local interstellar medium. Dosage time 

901 profiles are computed for isotropically incident fluxes onto flat surfaces from the flux spectra 

902 in Figure 3. The GEANT radiation transport code (http://wwwasd.web.cem.ch/ 

903 wwwasd/geant/), as implemented in our earlier work [Stumer et al., 2003; Cooper et al., 

904 2006], is used here for complete interactions of primary electrons, protons, and secondary 
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905 interaction products at energies above 10 keV from the flux spectra in Figure 3. Due to 

906 limitations on spatial step resolution in GEANT, the proton profiles have been extended to 

907 lower energies with stopping range and differential energy loss data of the Stopping and 

908 Range of Ions in Matter (SRIM) model [ Ziegler et al, 1985] (http://www.srim.org/) also used 

909 in the earlier work [Cooper et al., 2001]. 
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Figure 1 . Illustration of radiolytic gas-driven cryovolvanism for potential subsurface liquid 
water reservoir near the surface of Enceladus. A continuous rain of energetic electrons 
(yellow arrows) drives radiolysis and saturates the upper ice surface with oxidants (blue), 
mostly H2O2 but with mix of other species. Transient plume activity ejects ice grains falling 
back to the surface (white ballistic curves) and depositing multiple layers of water frost 
(white) interspersed with oxidant concentrations accumulating during lower plume activity. 
Ice upwelling continuously supplies fresh “fuels” (red), e.g. hydrates of CH4 and NH4, from 
the deep interior. Exposure to increasing temperature within the ice margins of the fluid 
reservoir initiates sequences of exothermic reactions from H2O2 dissociation to fuel 
oxidation, resultant gas production, and fluid reservoir heating. Percolation of expanding hot 
gas bubbles from the margins into the incompressible water fluid becomes the driving force 
for upward movement (white arrows) of gas-saturated fluid to form the plumes. Interactions 
of lower (left side) or higher (right side) concentrations of oxidants and fuels produce 
correspondingly less or more gas. 
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927 Figure 2. Representative differential flux spectra of Saturn magnetospheric electrons at the 

928 equatorial jovicentric orbits of the icy moons Mimas, Enceladus, Tethys, Dione, and Rhea as 

929 taken directly from the SATRAD model [Garrett et al., 2005], 
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Figure 3. Differential flux spectra of energetic electrons (red curves) at Enceladus and 
Europa, respectively from the planetary magnetospheres of Saturn and Jupiter, and of cosmic 
ray protons (black curves) irradiating icy bodies within and beyond the heliosphere. The 
thickest solid Enceladus curve is from the SATRAD model. The two medium-thickness 
Enceladus curves are respectively upper and lower limits from Voyager and Cassini/MIMI 
electron data. The upper limit comes from Voyager measurements compiled by Maurice et 
al. [1996], and the lower limit is from Cassini electron flux measurements by the 
Magnetospheric Imaging Instrument (MIMI) experiment [Krimigis et al., 2004] near the 
Jupiter-fixed longitude of this moon during the 9 March 2005 flyby. The limiting flux curves 
for Enceladus are extrapolated into the SATRAD curve at higher energies. The thinnest red 
curves are MIMI measurements for other selected periods of Enceladus orbit crossings: Julian 
days 48, 195, and 266 of 2005. A small peak at 0.03 MeV in the spectrum from day 48 is 
from a transient injection event. The Europa flux spectrum [Cooper et al., 2001] is a 
composite of Galileo Orbiter data below 1 MeV from the Energetic Particle Detector 
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946 [Williams et al., 1992] and Pioneer-Voyager model data [Divine and Garrett, 1983] at higher 

947 energies up to 40 MeV. The Enceladus and Europa model spectra are extrapolated for surface 

948 interaction modelling as power laws to 1000 MeV from the lower energy measurements. 

949 Proton flux spectra are also shown as derived for cosmic ray protons (four dashed curves) in 

950 the outer supersonic heliosphere within the classical Kuiper Belt near 40 AU and in the outer 

951 heliospheric boundary regions [Cooper et al., 2006] of the heliosheath and the local 

952 interstellar medium. The two 40-AU spectra are respectively for minimum [Cooper et al, 

953 2003] and maximum solar activity, the latter data being provided from previously 

954 unpublished Voyager data of co-author M. E. Hill. 
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Figure 4. Surface depth profiles in H 2 O ice for time in years to accumulate chemically 
significant dosages of 60 gigarad (~ 110 eV per irradiated H 2 O molecule) from total 
magnetospheric electron irradiation of Enceladus and Europa as compared to trans-neptunian 
objects. The two thin solid curves show the partial contributions of higher energy electrons at 
1-10 MeV and 10-100 MeV to the Enceladus total irradiation at 10 keV to 1000 MeV, as 
compared to the total irradiation profile for Europa (dot-dot-dash curve). Comparative 
profiles (dashed curves) are also shown for cosmic ray proton irradiation of objects near 40 
AU near solar minimum [ Cooper et al., 2003] and on highly eccentric orbits [ Cooper et al., 
2006] passing through the heliosheath and into the local interstellar medium. Dosage time 
profiles are computed for isotropically incident fluxes onto flat surfaces from the flux spectra 
in Figure 3. The GEANT radiation transport code (http://wwwasd.web.cem.ch/ 
wwwasd/geant/), as implemented in our earlier work [Stumer et al., 2003; Cooper et al.. 
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970 2006], is used here for complete interactions of primary electrons, protons, and secondary 

971 interaction products at energies above 10 keV from the flux spectra in Figure 3. Due to 

972 limitations on spatial step resolution in GEANT, the proton profiles have been extended to 

973 lower energies with stopping range and differential energy loss data of the Stopping and 

974 Range of Ions in Matter (SRIM) model [Ziegler et al., 1985] (http://www.srim.org/) also used 

97 5 in the earlier work [ Cooper et al . , 200 1 ] . 
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